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Abstract 

Using renormalization group with included self-energy effects and exact diagonalization of small 
clusters we investigate the ground state phase diagram of a two-dimensional extended Hubbard 
model with nearest-neighbor exchange interaction J, in addition to the local Coulomb repulsion U. 
The main instabilities are antiferromagnetism close to half- filling and d x 2_ y 2 superconductivity in 
the doped system. We find that self-energy effects are fatal for superconductivity in the repulsive 
Hubbard model (i.e. J = 0, U > 0). The superconductivity is triggered by finite J. The combined 
action of J and U interactions provide a remarkably efficient mechanism to enhance both d x 2_ y 2 
superconducting and antiferromagnetic correlations. 
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INTRODUCTION 



One of the most striking features observed in the phase diagram of the high T c super- 
conducting cuprates is the proximity between the insulating state with long-range anti- 
ferromagnetic order and the superconducting phase. This remarkable issue and the fact 
that at optimum doping the antiferromagnetic coherence length remains finite (£ ~ 4a for 
La2_a.Sr3.CuO4), have been the source of inspiration for many theoretical works in which the 
pairing mechanism is proposed to be originated in the fluctuations of the spin-density wave 
phase. Another peculiarity of these materials is that the superconducting order parameter 
has d x 2_ y 2 symmetry. The reason for this property is also thought to be the closeness to the 
antiferromagnetic phase as the lines of nodes of the superconducting gap allows for the ex- 
istence of gapless spin excitations and superconductivity can coexist with spin-density wave 
fluctuations. Very recently, the proximity between these two phases supplemented with the 
effect of disorder has been the basic ingredient in a phenomenological model to suggest that 
colossal effects can be expected in the phase diagram of the cuprates 

The fact that the ground-state of the undoped materials is antiferromagnetic immediately 
suggests that the Hubbard model could be a good candidate for a microscopic description 
of these compounds. Instead, the explanation of the superconducting mechanism in the 
framework of this model remains very controversial. A superconducting solution is found 
when the model is tackled with some many-body techniques (2), U, y, _3, Q, _1 _1 __|. 
However, all the numerical works devoted to search indications of superconductivity in the 
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HQ, Q. 



Another candidate 



repulsive Hubbard model have been negative so far 
model to provide the basis for the theoretical investigation of the high T c superconductors 
is the t — J model, which coincides with the Hubbard Hamiltonian in the limit of U — > 00 
(corresponding to J — > 0). In the case of the t — J model, many-body methods and numerical 
results seem to be in agreement concerning the possibility of a superconducting state with 
d x 2_ y 2 symmetry Q, Q, [3]. The range of J at which it may occur is, however, not 
precisely determined and it is likely to lie beyond the region where the mapping from the 
Hubbard model is valid. 

This motivates the study of the t — J — U model 

H = -t^2 ( c l* c jv + hc) + S i' S i + U J_ n i\ n Hi (!) 
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which, in addition to the Coulomb repulsion U of the usual Hubbard Hamiltonian, contains 
a nearest-neighbor exchange interaction J as the t — J model. In such a way, we can expect 
to retain appealing features of both models, like the charge fluctuations introduced by U but 
forbidden in the constrained t — J model and the robust superconducting correlations intro- 
duced by the exchange interaction and explore the interplay between both effects. Beside 
these heuristic arguments, the t — J — U model is closely related to the extended Hubbard 
model with correlated hopping, which has been derived from the three-band extended Hub- 
bard model as an effective one-band Hamiltonian to describe the low energy properties of 
the Cu-0 planes of the high-T c materials Q|. The latter model is characterized by three 
different nearest-neighbor hopping amplitudes t A Ai ^ab-, ^bb depending on the occupation of 
the two sites involved, as well as the usual Coulomb repulsion U, 

H = U^2 n^mi - ( c Lcja + he) 

{*aa(1 - n iw )(l -n jw ) 

t AB [(l - n iw )n jw ) + n iw (l - n jw )} 

t B Bn iw n jw }. (2) 

Each of the parameters of the above model depends on the parameters of the original three- 
band Hamiltonian and there is a degree of uncertainty in their precise values. The rele- 
vant property is that reasonable estimates indicate that £ab is larger than the other two 
|ig| | . The analysis of the different hopping processes in (J2J) reveals that the one driven 
by t A B mediates antiferromagnetic correlations. In particular, in the strong coupling limit 
U » tAA,tAB,tBB, the exchange interaction obtained by treating ((2)) with second order 
perturbation theory is J = 4t AB /U. Furthermore, when this process is suppressed, antifer- 
romagnetic correlations are completely eliminated and a metal-insulator transition can take 



place at finite U in the half-filled system [20, |2jJ, |22] . For weak coupling, the Hamiltonian 



ce techniques in 2D 



23j and 



(J2J) can be treated with mean-field Hartree-Fock and BCS-li 

with operator product expansion and bosonization in ID [24 1. It turns out that, close to 
half-filling, the predicted phase diagram is equivalent to that obtained starting from an ef- 
fective t — J — U model 25j] . In summary, the relevance of the t — J — U model to describe 
the Physics of the cuprates can be also supported by its closeness to the correlated Hubbard 
model derived from the more detailed three-band Hubbard model for the Cu-0 planes. 
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Coming back to the phenomenology of the superconducting cuprates, the t — J — U model 
provides the framework for recent suggestions based on a combined order parameter with 
superconducting and bond- density- wave components with d x 2_ y 2 symmetries to explain the 
intriguing pseudogap phenomena [2^. In general, extended Hubbard models in low dimen- 



sions are of interest in the context of high-T c materials 



27, 



28, 



and also to understand 



the rich structure observed in the phase diagrams of organic materials like (TMTSF^PFg, 



and (TMTSF) 2 C10 4 
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which can not be explained on 



the basis the local Coulomb interaction of the usual Hubbard model. 

The t — J — U model has an explicit effective attraction mediated by J that enables 
a superconducting solution even within a simple BCS-like description. Numerical studies 
of this model in ladders indicate that superconducting correlations with d x 2_ y 2 symmetry 
are enhanced in comparison to those of the usual Hubbard model Q|. More recently, an 
analysis based on the Landau Fermi liquid picture has been adopted to argue that the pairing 
interaction mediated by J combined by a strong renormalization of the effective density 
of states caused by U results in a significant enhancement of the superconducting order 
parameter while quantum Monte Carlo simulations based on a variational wave function 
with a BCS structure support this picture j^. Investigations of the weak coupling phase 
diagram at van Hove fillings of the two dimensional version of this model supplemented by 
a nearest-neighbor Coulomb repulsion and next-nearest neighbor hopping amplitude have 
also been reported j^jj]. 

The aim of this work is to investigate the two-dimensional phase diagram of the t — J — U 
model with two complementary techniques: exact diagonalization (ED) of a small cluster 
with 4x4 sites and the one-loop renormalization group (RG) technique presented in Refs. 
Q| • The latter are expected to provide reliable indications of the main instabilities of the 
Fermi liquid in the limit of weak interactions. ED is an unbiased method that consists in 
the exact calculation of the ground state wave function, allowing for the direct evaluation of 
the relevant correlation functions but has the drawback that the size of the clusters that are 
amenable to be treated is small. In spite of that limitation, in the case of the t — J model, 

results on such small clusters are among the most robust 



the conclusions based on numerica 
ones regarding superconductivity 17]. It is therefore interesting to analyze the predictions 
of these two methods and to compare them with previous results. 

The article is organized as follows, in section II we provide some technical details on the 
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methods we employ. In particular, we describe the procedure followed to include self energy 
corrections in the RG treatment. An important reference is the behavior of the relevant 
susceptibilities evaluated with this RG procedure for the usual Hubbard model. This is 
discussed in Section III. In sections IV and V we present the results for the t — J — U model 
using RG and ED, respectively. Section VI is devoted to summary and conclusions. 



TECHNICAL DETAILS 



Renormalization group method 

The basic hypothesis of the renormalization group method is that the normal state is well 
described by the effective action for quasiparticles near the Fermi surface. The description of 
the many body problem is done in terms of the two-body effective interaction Ui(9i, 62, #3) 
and the quasiparticle weight Zi{9). In the Wilson's renormalization scheme, £^(#1, #2, #3) 
is the effective interaction between electrons within the ring ±A around the Fermi surface. 
This interaction is renormalized by the scattering processes involving all electrons outside 
the ring ±A. A physical interpretation of the cutoff A is that it plays the role of an effective 
temperature or the experimental probe frequency. In previous versions of RG 0, 0, 0, , 
it was assumed Zi(6)=l. In what follows, we summarize the improved RG method of Ref. 



41 1 which also considers self-energy corrections by taking into account the renormalization 
of the quasiparticle weight. 

An important issue to note is the fact that 1/1(81,62,63) and Z\{&) depend only on the 
angles &i that parametrize the positions of the particles on the Fermi surface. This is 
justified by a simple power counting which tells us that only the angular dependence of 
the effective interaction is marginal (or marginally-relevant) and only terms up to linear 



in energy are to be kept in the renormalization of the angle-dependent self-energy 41]. 
In the renormalization group procedure, these functions are continuously renormalized as 
the energy cutoff, parametrized by the scale I as A = 8texp(— I) is reduced. The ensuing 
equation for the evolution of Ui within the one-loop renormalization group scheme has the 
following structure: 
BU 

-Q+ = Ppp{U h U t } + 2(3 ph {U h U t } - /3 ph {U h XUi} - f3 ph {XU h U t } - Xf3 ph {XU h XU t }, (3) 
where X is the exchange operator defined as XU(1, 2, 3) = U(2, 1, 3). One must remember 




FIG. 1: From the left: Particle-particle (Cooper) and particle-hole (density-wave) differential 
bubbles (3 pp (0i, 02, #3) and (3 p h(0i, O2, #3). Propagators corresponding to degrees of freedom at 
the Fermi surface and outside the ring ±A are indicated with a bar and a >, respectively. 

that © is a functional flow equation, i.e. Ui and all terms on the right-hand side depend 
on three angles (9±, 9 2 , 9 3 ). Particle-particle (Cooper) and particle-hole (density-wave) dif- 
ferential bubbles j3 pp and (3 p h are shown on figHJ To solve numerically Eq.(j3J) we discretize 
the angle 9 which defines the so-called iV-patch model. The function Ui(6i, 92, #3) is then 
represented by a set of coupling constants labeled by three discrete indexes. For the t — J — U 
model, the initial condition is 

Ui =0 (9!, 9 2 , 9 3 ) = |[/ - ^[cosfo - + cos(y 3 - y x )) - ^[cos(x 3 - x 2 ) + cos(y 3 - y 2 )\ 

'(4) 

where (xi,yi) are components of the wave vector kj, corresponding to the angle 9\ at the 
noninteracting Fermi surface. All coupling constants are found to diverge at the same critical 
scale l c like 

where the weights U are model dependent constants. This type of solution is called the 
fixed-pole solution in contrast to the mobile-pole solution, where different coupling con- 
stants diverge at different critical scales. For realistic systems, where the initial coupling is 
not extremely small, only fixed poles are relevant [4^. The critical scale l c depends on the 
bare coupling constants U and J and on the band filling parametrized by the chemical poten- 
tial. The critical cutoff A c = 8texp (— l c ) appears to be the characteristic temperature of the 
model. The most precise non-restrictive interpretation of A c is that at this energy electrons 
start to build bound states. Namely, the poles in Eq|S]are the two-particle propagator poles, 
indicating the on-set of bound states. These bound states can be of charge, spin or super- 



conducting kind and all of them renormalize the one-particle weight. This renormalization 
is angle dependent and its evolution is described by the following expression 

dtiogztf) = J de' j{e',-K) m {e,e') , (6) 

with the initial condition Z\{9) = 1. The function J7(0, e) is the angle dependent density 
of states at the energy e (measured from the Fermi level). The quantity rji(9,9') contains 
particle-particle (pp) and particle-hole (ph) contributions 

7j,(0, 9') = {2X-l)f3 pp {U h U l }+2(3 ph {XU l , XU l }+2f3 ph {U h *7,}-/V{^, XU l }-f3 ph {XU u 17,} 

(7) 

with all terms on the right-hand side taken with external legs (01,02,03) = (9, 9', 9'). The 
interaction inserted in all beta functions obeys the scaling equation (JBJ). The relation between 



interaction Ui and the usual physical interaction F contains the rescaling Z-factors, 



4l|: 



r,(l, 2, 3) = [Z,(1)Z,(2)Z,(3)Z,(4)]- 1 / 2 C/,(1, 2, 3) . (8) 

To find out which correlations are relevant and candidates for order parameter we must 
allow the theory to choose between all possible 2-particle correlations. For this reason 
we have to follow the renormalization of several angle-resolved correlation functions. The 
superconducting correlation function xf C {@i>@2) measures correlations between the Cooper 
pairs (9\, 0i + 7r) and (6*2, 6*2 + 71"), all states being at the Fermi surface. The antiferromagnetic 
correlation function x? F (@i>@2) correlates two nested electron-hole pairs c Le 1 ) crc k(0i)+(7r,7r) 
and c^g^o" Ck(e 2 )+(7r,7r)- The charge density wave correlation function x? DW i.^1^2) correlates 
the nested charge-like electron-hole pairs. To get the renormalization group flow of all 
correlation functions we follow the procedure given in Ref.jfJ but dressing the electronic 
propagators with Z-factors as in Ref.j^lj]. We get 

Xf (01, 62) = Zl{ QJ Zl{ Q 2) f ^ ~ Z l^ 9 ) A 5 (0)^(0, *2) • (9) 

The function Df(9) writes 

Df c {e) = \ J2 J&HW) (io) 



for the SC channel and 



^ U(-A(i),fl) 
Dl W "21 + H/A(/)' 



for the AF channel where only the negative shell (y = —1) contributes to the flow. 
The flow of the quantity zf(0i, 6) that has the role of a triangular vertex writes: 

[di - v(0i) - v(o 2 )] o 2 ) = -j do z*{o lt e)D & l {6)v l 5 {e, e 2 ) . (12) 

The meaning of zf(0i, 9 2 ) is that 

zf (6 U 6 2 ) = Z l (0 1 )zf(0 1 ,0 2 )Z l (0 2 ) 
so that the initial conditions for zf and for zf are the same: 

zf =0 (e 1 ,e 2 ) = 5 D (e 1 -e 2 ), (13) 

where 5d is the Dirac function. Initial conditions for susceptibilities are: 

xlo(^iA) = 0. (14) 

After discretization we integrate numerically equations and (|12p. The relevant suscepti- 
bility in each channel is the dominant eigenvalue of the angle-resolved correlation function. 
The corresponding eigenvector determines the angular dependence of the order parameter. 

Exact diagonalization 

This method consists in the exact computation, by recourse to Lanczos algorithm, of 
the ground state (GS) wave function of the model Hamiltonian (0) on a small cluster. We 
consider a cluster containing 4x4 lattice sites. 

Due to the finite size of the cluster, its reciprocal lattice contains just some few (N = 16) 
points. As several instabilities are expected to be competitive in this model, the rough mesh 
of available k-vectors can produce an important bias and lead to an artificial enhancement of 
some kind of correlation. The only resort to minimize this drawback is to consider different 
kinds of boundary conditions, which is equivalent to consider different choices of the 16 
k-points. An arbitrary choice of boundary conditions breaks the point group symmetry 
of the original lattice. Exceptions are the periodic (P) and antiperiodic (AP) boundary 
conditions, which lead to the mesh indicated in FigsOt 1 and c, respectively. Mixed (M) 
boundary conditions (periodic in one direction and antiperiodic in the other) lead to the 
pattern depicted in FigJ2^. The latter breaks some of the symmetries of the C^ v group of 
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FIG. 2: Reciprocal lattice of the 4x4 cluster assuming (a) mixed boundary conditions, (b) peri- 
odic boundary conditions and (c) antiperiodic boundary conditions. Ground state configurations 
corresponding to the fillings n = 0.75 (thick arrows) and n = 0.875 (adding the thin arrows to the 
previous case) are also sketched. 

the original lattice, the corresponding point group being Civ- Typically, these three choices 
of boundary conditions are the ones leading to the lowest energy. In our study, we compute 
the GS of (|T|). considering all the three possibilities above mentioned within the subspaces 
corresponding to the different one-dimensional representations of the point group and total 
k = (0,0). In the noninteracting system, it is easy to see that closed-shell configurations 
are those leading to the lowest energy. For some densities of particles this condition is, 
however, not possible to be fulfilled for any choice of the boundary conditions and the GS 
is degenerate. The interactions normally lift most of the degeneracies. In some cases, it is 
observed that when the interactions overcome some particular value, a change is produced 
in the BC leading to the lowest energy. The latter effect is an indication that interactions 
lead to some qualitative change in the behavior of the GS. This is, of course, a mere finite 
size effect but provides a valuable information, since it reflects that the system prefers a 
change in the population of available k-points in order to take advantage of the interactions 
and thus lower the energy. 

To investigate the superconducting correlations in the GS it is useful to study the behavior 
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of the pair correlation function (PCF) 



P(r) 



1 



^(^olAKRi + rJA^Ri)!^), 



(15) 



N 



where |\I/o) is the ground-state wave function while A£ s (Rj) = cLcj, for on-site s pairing, 
A^(Rj) = J2s f^)[4+S]4i ~ 4+54^]/^' witn fes(5) = 1 for extended s pairing, and 
fd(8) = 1 (fd($) — — 1 ) when 5 = ±(1,0) (5 = ±(0,1)) for d x 2_ y 2 pairing. This function 
is normalized in such a way that |A^(i)|\l/o)| 2 = 1. A superconducting state with pairs of 
a given symmetry is expected to have sizable correlations between pairs far separated by 
arbitrary large distances. In the case of the 4x4 cluster, the largest available distance is 
r = (2,2). The PCF between pairs separated by this maximum distance is denoted P m . 
As even in the noninteracting limit the PCF can be finite, we interpret an enhancement 
of the corresponding PCF relative to its value at U = 0, J = as an indication of the 
superconducting instability. 

To study the spin-density-wave (SDW) correlations it is useful to compute the spin-spin 
correlation function, 



which provides information on the nature of the correlations between spins. 

RENORMALIZATION-GROUP RESULTS FOR HUBBARD MODEL WITH RE- 
PULSIVE AND ATTRACTIVE INTERACTION 

In this section we present results for the usual Hubbard model obtained by the RG method 
described in the previous section. The motivation is twofold. For the case of the U > 
model, it was shown in Ref. j^j that self-energy corrections included in the renornalization 
of Z(9) are important to predict the antiferro magnetic instability at half-filling in 2D and 
to recover the correct expression for the jump at the Fermi points in ID but the behavior 
of the susceptibilities away from half-filling has not been analyzed so far. On another hand, 




(16) 



and to analyze the Fourier transform 




(17) 
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FIG. 3: (Color online) AF (upper panels) and d x 2_ y 2 SC (lower panels) susceptibilities as functions 
of the scale I for Hubbard model with U = 1.6i and densities determined by = oo (left-hand 
panels) and = 3 (right-hand panels), corresponding, respectively, to half-filling (parquet regime) 
and high doping (within the BCS regime). Thick and thin lines correspond to results with and 
without self-energy corrections, respectively. The noninter acting particle-particle and particle-hole 
susceptibilities are shown in dashed lines. The critical scale l c as well as the scale I* at which the 
susceptibility departs from the behavior of the noninteracting one are indicated with magenta and 
cyan arrows, respectively. 

for U < the model is a paradigmatic example of a superconductor and it is, therefore, an 
important reference point to analyze the behavior of the superconducting correlations. The 
results shown correspond to a discretization of 32 patches. 

When applied to the usual Hubbard model with repulsive interactions (U > 0), the RG 



without considering self-energy corrections 



predicts that at some critical scale l c , 



the superconducting susceptibility with d x 2_ y 2 symmetry together with AF susceptibility 
diverge within the parquet regime < A), even at half- filling. This feature is expected 
to be an artifact of the approximation since for n = 1 umklapp processes are active and 



they are expected to drive the system toward an insulating AF state. In Ref. |4l| it was 
shown that, due to the renormalization of the quasiparticle weight, the SC susceptibility is 
suppressed below its noninteracting value, while the AF susceptibility remains enhanced, 
although weakly, relative to the noninteracting particle-hole one. 

Fig. El shows the behavior for the different susceptibilities as functions of the scale /, at 
different densities denned by the chemical potential \i = 8texp(— Z M ). The susceptibilities 
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with and without self-energy corrections are respectively plotted in thin and thick lines. 
The two left-hand panels correspond to half-filling (7 M = oo) and summarize the results of 
Ref. j^lf. Without self-energy corrections both SC and AF susceptibilities diverge at l c , AF 
being dominant. Self-energy effects suppress the divergence, decreasing the SC susceptibility 
bellow its noninteracting value. The AF susceptibility is also renormalized by self-energy 
effects but remains larger than the noninteracting one. The scale I* (see Fig. |3J) at which the 
AF susceptibility begins to depart from the behavior of the noninteracting one is, however, 
not affected by the self-energy corrections. This kind of behavior extends along the parquet 
regime, defined by the condition |/i| < A. 

The region at lower densities such that > A is usually called BCS regime because 
only Cooper channel has logarithmic contributions to the effective interaction. The two 
right-hand panels of Fig. El show how in this regime both AF and SC susceptibilities remain 
slightly larger than their noninteracting values when self-energy corrections are taken into 
account. 

These results need some discussion. Up to now it was rather widely accepted that the 
2D Hubbard model has a d x 2_ y 2 SC state away from half-filling. On the contrary results 
on the right-hand panels of Fig. |3] do not convincingly indicate a strong superconducting 
instability. The SC correlation functions being only weakly enhanced and the scale I* being 
strongly renormalized towards high values when self-energy effects are taken into account 
(see right-hand lower panel of Fig. |3J). This indicates that the energy A* = 8texp(— /*) 
at which superconducting correlations begin to manifest themselves decreases when self- 
energy effects are considered. At A c there is a creation of bound states because the 4-points 
vertex (j3J) has poles. These bound states are the gapless modes that destroy the one-particle 
coherence via self-energy corrections and the result is that the phase transition is suppressed, 
probably to some finite lower energy. 

To justify the above interpretation we calculated for comparison the RG flow of the at- 
tractive Hubbard model, for which the instabilities are better known. The U < model 
is a paradigmatic example of s-wave type superconductor and its phase diagram has been 
investigated in detail by mean-field and numerical techniques. Unlike the repulsive coun- 
terpart, predictions by different methods agree about the main instabilities expected in its 



phase diagram 



44j : At weak coupling, the explicit local negative interaction leads to 



BCS-like superconductivity away from half-filling. For n = 1, local CDW is believed to be 
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FIG. 4: (Color online) CDW (upper panels) and s-wave SC (lower panels) susceptibilities as 
functions of the scale I for U = — 0.8t and densities determined by = oo,5,3, corresponding to 
half-filling and two densities within the BCS regime (from left to right). Other details are the same 
as in Fig. [21 

I — I 

degenerate with s-wave superconductivity in the GS [43]. In ID, such degeneracy is exact 
due to symmetry reasons. This is because the repulsive model has dominant AF correlations 
with a power law decay. Since no breaking of the SU(2) symmetry is possible, the behavior 
of the correlations in any of the spacial direction must behave in the same way. On another 
hand, there is an exact transformation cL — > (— l) l Cjf, c|, — > ct, which maps the repulsive 
model into the attractive one, while maps the degenerate z and x, y components of the 
dominant SDW correlations of the U > model to the CDW and superconducting ones, 
respectively, of the U < counterpart. The above reasoning can also be extended to the 
2D case provided that no symmetry breaking in the ground state of the U > model takes 
place. Results obtained with RG for the attractive model in 2D shown in Fig 0] are in very 



43, 



4J. At half-filling 



good agreement with the description provided by numerical methods 
(In = oo, see left-hand panels of Fig. HJ), the most remarkable feature is the degeneracy 
observed between CDW and SC susceptibility with local s-wave symmetry, which remain 
slightly larger than the noninteracting one when self-energy corrections are considered. This 
should be rather expected since the present method provides a description of the normal 
state and only the on-set of the instability towards the symmetry broken state is captured. 
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Bellow half-filling (Z M = 5, 3, cf. middle and right-hand panels of FigHJ), CDW susceptibility 
becomes weaker, approximately equal to the noninteracting one, while the superconducting 
susceptibility becomes more enhanced. Self-energy effects suppress the divergence of the 
different instabilities but do not renormalize the scale I*. 

Bellow half-filling, the s-wave type SC susceptibility is clearly larger than the noninter- 
acting particle-particle one. As expected, the CDW remains non renormalized. Compared 
to the repulsive Hubbard model (Fig|3J), the enhancement of SC correlations is more con- 
vincing and in the present case we can indeed assign the divergence at l c to the on-set of the 
s-wave superconductivity. Altogether, the attractive Hubbard model away from half-filling 
shows strong tendencies towards the s-wave superconductivity even when the self-energy 
corrections are taken into account. This is in contrast to the repulsive case where the self- 
energy effects have a more pronounced effect against the d x 2_ y 2 superconductivity. The 
reason for this behavior is rather simple: the s-wave superconductivity in the attractive 
Hubbard model is an effect of first order in U, while the self-energy corrections are of second 
order. As long as the coupling is weak the BCS like instability is a good approximation. 
Formally this means that the mean-field and the Kosterlitz-Thouless transitions are close to 
one another. In the case of the repulsive interaction the d x 2_ y 2 superconductivity and the 
self-energy corrections are both of second order in U, so that by decreasing U one cannot 
promote only superconductivity and make the fluctuations subdominant. 

RENORMALIZATION-GROUP RESULTS FOR THE t- J-U MODEL. 

We present results for the relevant susceptibilities in the t — J — U model in Fig. [SJ 
The left-hand panels correspond to half-filling and the behavior is similar to that of the 
Hubbard case. Namely, SC susceptibility with d x 2_ y 2 symmetry evolves to values below of 
the noninteracting one, while AF susceptibility remains higher than its noninteracting value, 
indicating that at n — 1, the system flows towards an AF insulating state. Upon doping the 
d x 2_ y 2 superconductivity gets progressively stronger and becomes the dominant instability 
at about the crossover line A = fi, just as in the Hubbard model. The right-hand panels 
show a typical flow of the susceptibility in the BCS regime (A < /i) where the dominant 
correlations are of superconductivity type. 

The above picture shares some features with the behavior observed in the Hubbard model, 
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FIG. 5: (Color online) AF (upper panels) and d x 2_ y 2 SC (lower panels) susceptibilities as functions 
of the scale / for Hubbard model with U = J = 1.6i and densities determined by = oo (left-hand 
panels) and lp = 3 (right-hand panels), corresponding, respectively, to half-filling (parquet regime) 
and high doping (within the BCS regime). Other details are the same as in Fig. |31 

discussed in the previous section. Important issues to highlight are: (i) At half-filling the 
behavior of the relevant susceptibilities is very similar to that of the repulsive Hubbard 
model. However, the scale l c at which the on-set of the AF instability takes place as well 
as the scale I* are smaller than for the J = case. The susceptibility is also significantly 
larger than the noninteracting one, even when self-energy corrections are included in the 
RG procedure. These features indicate that J contributes to increase the AF correlations 
and the Neel temperature, (ii) At higher dopings, within the BCS regime, AF susceptibility 
coincides with the noninteracting one while SC correlations become significantly enhanced. 
This is in contrast to the behavior of the repulsive model (cf. Fig. EJ) where SC correlations 
are only weakly enhanced. Instead, the behavior on the right panels of Fig. 03 resembles the 
one of the attractive Hubbard model, if we associate d x i_ y i SC and AF susceptibilities in 
Fig. |3] respectively to s-wave SC and CDW susceptibilities on middle and right panels of Fig. 
0] Also note that, as in the attractive model, the scale I* remains unaffected by self-energy 
effects. In addition, the scale l c at which the on-set of the superconducting instability is 
observed is small in the t — J — U model, implying a high critical temperature. We have 
carried out a similar analysis for other values of the parameters J and U and found that the 
symmetry of the dominant superconducting correlations is always d x 2_ y 2. We have verified 
the reliability of these results upon increasing number of patches up to 64 patches. The 
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FIG. 6: (Color online) Phase diagram showing combined action of U and J interactions. The 
"critical temperature" is defined as T c ~ A c = 8texp(— l c ) (see text). Left and right panels 
correspond to J = 1.6t and J = OAt, respectively. Plots in blue circles correspond to U = 
0, 0.4i, 0.8t, t, 1.6t, 3.2t. The ones corresponding to the lowest and highest U are drawn in open 
and dark filled symbols, respectively. The plot in red squares corresponds to the usual Hubbard 
model (J = and U = 1.6t). The line T = (J, separates the regions where the SDW and d x 2_ y 2 
superconducting (dSC) instabilities are the dominant ones. 

reason for the robustness of the SC correlations is in the fact that the J interaction has an 
attractive d x 2_ y 2 SC component, so that superconductivity exists already at first order of U\ 
while the fluctuations are only subdominant, just as in the attractive Hubbard case. 

In order to have a more quantitative representation of the role played by both interactions 
we present the phase diagram in Fig. El The "critical temperature" plotted in the figure is 
defined as T c ~ A c where A c = 8texp(— l c ). Results for the repulsive Hubbard model are 
also shown for comparison. By comparing the plots for J = 0, U = l.Qt and J = 1.6t, U = 
in the left panel and for J = 0, U = 1.6t and J = 0At,U = 1.6t in the right panel (the 
arrows are drawn to ease the reading), it is clear that J is a remarkably efficient mechanism 
to drive AF close to half- filling and d x 2_ y 2 superconductivity in the doped system. Another 
important feature is that at fixed J, the effect of U is to increase the "critical temperature" . 
This means that the two interactions are not competitive but, instead, cooperate to increase 
the strength of antiferromagnetic and superconducting correlations. 
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EXACT-DIAGONALIZATION RESULTS FOR THE t-J-U MODEL IN THE 4x4 
CLUSTER. 

The results of the previous section suggest that the combined effect of the interactions 
J and U drives superconductivity with d x 2_ y 2 symmetry in the doped system, leading to a 
significant enhancement of the superconducting correlations. As in the case of the pure Hub- 
bard model numerical methods fail to detect the tendency towards d^^-superconductivity, 
it is interesting to analyze the case of finite J. 

We show bellow results obtained by following the strategy explained in Section II. We 
computed the exact GS energy and wave function in a 4 x 4 cluster and calculated the 
correlation functions between pairs with local and extended s-wave and d x 2_ y 2 symmetries. 
We focused our attention on the study of three different fillings (n = 0.625,0.75,0.875), 
corresponding to N = 10, 12, 14 particles in the cluster. 

Let us begin with the analysis of the boundary conditions leading to the lowest energy. For 
the case n = 0.625, the GS is obtained for PBC. In the noninteracting limit, it corresponds 
to a closed-shell configuration and for all the explored values of the interactions, it lies in 
the subspace associated to the representation of the point group with s-wave like character. 
The k-points lying on the Fermi surface of the noninteracting system are (tt/2,0) and the 
symmetry related points. For these points, the structure factors / es (k) = cos(k x ) + cos(k y ) 
and /d(k) = cos(k x ) — cos(k y ), corresponding to BCS gaps with extended s and d^.^-wave 
symmetries have exactly the same strength |/^| = \f es \ = 1. 

For fillings n = 0.75 and n = 0.875, the lowest energy in the noninteracting case is 
achieved by considering MBC. In this limit, the only ingredient playing a role in the energetic 
balance is the kinetic energy gain. As interactions are switched on, the GS corresponds to 
APBC for sufficiently large J and U. This indicates that the k-points tuned by APBC are 
able to take advantage of some effect of the interactions, compensating the loss of kinetic 
energy. For the latter boundary conditions, the Fermi points of the noninteracting system 
lie on the lines of nodes of f es . Therefore, it is likely that the most favored instability by 
such a change of population in the k-space is d x 2_ y 2-wave superconductivity. Some values 
of the GS energy per site are shown in tables HI and ITT1 

It is interesting to note that, for the density n = 0.75, the GS belongs to the representation 
of the point symmetry group with s-wave-like character within the region where the GS 
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.7092 
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.3651 
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.4584 


-1 


.5100 
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-1.6864 


-1.7526 


4 


-1. 


1607 


-1 
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-1 


.2616 


-1 


.3176 


-1 


.3772 


-1.5075 


-1.5803* 


-1.6619* 
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-1. 
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-1 


.1045 


-1 


.1615 


-1 


.2225 


-1 


.2874 


-1.4291* 


-1.5125* 


-1.5975* 
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-0. 


.9774 


-1 


.0326 


-1 


.0924 


-1 


.1565 


-1 


.2247 


-1.3773* 


-1.4630* 


-1.5503* 
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.9255 
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-1.3385* 


-1.4258* 


-1.5145* 
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.0735 
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-1.3087* 


-1.3969* 


-1.4867* 


14 
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.9786 
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.1187 


-1.2852* 


-1.3741* 
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.8377 
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.9567 
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-1.2662* 


-1.3556* 


-1.4465* 
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.8204 
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.9391 


-1 


.0073 


-1 


.0807 


-1.2505* 


-1.3403* 


-1.4315* 
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-0. 


.8063 
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.8621 


-0 


.9246 


-0 


.9929 


-1 
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-1.2375* 


-1.3275* 


-1.4189* 



TABLE I: Ground state energy per site for the 4x4 cluster with particle density n = 0.75. Stars 
indicate states with APBC and the representation of the point group that transforms like s-wave. 
Otherwise the states correspond to MBC. 

corresponds to APBC. Instead, for n = 0.875 and also within the region of parameters where 
the GS corresponds to APBC, the character of the point-group representation is d x 2_ y 2-wave 
like. Since in the present cluster n = 0.875, differs from n = 0.75 in two particles, this 
change of representation is consistent with the idea that a pair of particles with d x 2_ y 2- 
wave symmetry was added to a many-particle background with total s-wave symmetry. We 
actually speculate that such a background is also made up of paired particles. 

The behavior of P m , the pair correlation function ()15|) corresponding to pairs separated 
by the maximum possible distance of the cluster, for a particle density n = 0.625 is shown 
in FigEI The correlation function corresponding to local pairs with s-wave symmetry is 
much weaker and is not shown. The corresponding values for the noninteracting system 
are indicated in dot-dashed lines to ease the comparison. In the latter limit, correlations 
of pairs with d x 2_ y 2 symmetry remain weaker than those of the noninteracting case while 
s-wave ones are slightly enhanced for small enough U. The effect of J is to produce a weak 
enhancement of P m within the two symmetry channels in comparison to the pure Hubbard 
case. In particular, d x 2_ y 2 ones become stronger than those of the noninteracting case for 
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-1.0445 
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-0.7712 


-0.8580 


-0.9540 


-1.0566 


-1.1685* 


10 


-0.6978 


-0.7897 


-0.8908 


-0.9991* 


-1.1154* 


12 


-0.6460 


-0.7406 


-0.8446 


-0.9568* 


-1.0751* 


14 


-0.6079 


-0.7037 


-0.8095 


-0.9239* 


-1.0437* 


16 


-0.5788 


-0.6752 


-0.7819 


-0.8979* 


-1.0186* 


18 


-0.5560 


-0.6524 


-0.7598 


-0.8768* 


-0.9982* 


20 


-0.5376 


-0.6338 


-0.7417 


-0.8594* 


-0.9812* 



TABLE II: Ground state energy per site for the 4x4 cluster with particle density n = 0.875. 
Stars indicate states with APBC and the representation of the point group that transforms like 
d x 2_ y 2-wave. Otherwise the states correspond to MBC. 
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FIG. 7: (Color online) Pair correlation function at the maximum distance, P m , of the cluster for 
a density of particles n = 0.625. Blue circles correspond to pairs with d x 2_ y 2 and red diamonds 
to extended s symmetry. Different plots correspond to J = 0,0.25,0.5,0.75,1. Open and filled 
dark symbols correspond to the lowest and highest value of J, respectively. The dotted-dashed line 
indicates the value of P m in the noninteracting limit. 
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FIG. 8: (Color online) Pair correlation function at the maximum distance of the 4x4 cluster, P m , 
for a density of particles n = 0.75. Different plots correspond to J = 0, 0.25, 0.5, 0.75, 1, 1.5, 1.75, 2. 
Upper and lower panels correspond to APB and MBC respectively. Other details are as in Fig. [Tj 

small enough U. 

The pairing function P m for the density n = 0.75 is shown in Fig|H| We have analyzed 
the behavior using the two boundary conditions leading to the lowest energy. In the case of 
APBC shown in the upper panel of Fig. |SJ only the correlation of pairs with d x 2_ y 2 symmetry 
is shown, since those with local s- and extended s-wave are negligibly small in comparison. 
In the case of MBC shown in the lower panel, correlations in both, extended s-wave and 
c? x 2_ y 2-symmetry channels are only slightly enhanced for some values of J and suppressed 
for others. Instead, for APBC, a clear enhancement of the correlations of pairs with d x 2_ y 2 
is observed as J is switched on. Note that, in contrast to the case J = 0, the correlation 
function P m lies above the line indicating the magnitude of P m in the noninteracting case. 

Similar remarks apply to the behavior of P m at the density n = 0.875, shown in FigJUJ 
The enhancement of d x 2_ y 2 pairing correlations for the case of APBC is related to the fact 
that the available k vectors mainly populated in the ground state contribute with a sizable 
structure factor to a pairing interaction at the Fermi surface with d x 2_ y 2 symmetry. 

If, for the latter densities, we plot the pair correlation function in the GS corresponding 
to the optimal boundary condition (i.e. that leading to the lowest energy), we obtain the 
picture shown in Fig. in good agreement with the analysis done in the discussion about 
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FIG. 9: (Color online) Pair correlation function at the maximum distance of the 4x4 cluster, P m , 
for a density of particles n = 0.875. Other details are as in Fig. El 
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FIG. 10: (Color online) Pair correlation function at the maximum distance of the 4x4 cluster, 
P m , for pairs with d x 2_ y i symmetry corresponding to the optimal boundary conditions for the 
ground state. Upper and lower panels correspond to densities of particles n = 0.75 and n = 0.875, 
respectively. Other details are as in Figs. 0and|HJ 
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FIG. 11: (Color online) Fourier transform of the spin-spin correlation function for n = 0.75 and 
MBC. Different panels correspond to different U. Open and filled dark symbols correspond to 
J = and J = 1, respectively. Other plots correspond to intermediate equally spaced values of J. 
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FIG. 12: (Color online) The same as Fig. HI for n = 0.875. 



the behavior of the GS energy, we see that the change in the boundary condition leading to 
the lowest energy, is accompanied with an enhancement of the d^^-wave pairing correlation 
function. The pairing correlation functions with extended s symmetry are, instead, vanishing 
small within all the range of parameters. This could be an unfortunate consequence of the 
small size of the cluster and to the fact that the most populated k-points when the boundary 
conditions change, lie on the lines of nodes of the structure factor f es . However, this behavior 
is in agreement with the results predicted by RG. 

To finalize, we present some results on the behavior of spin-spin correlation functions. 
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Fig. and Fig ^] show some typical plots corresponding to the GS in the cluster with 
MBC. In the case of Fig. [TP the latter boundary condition corresponds to the one leading 
to the lowest energy within the range of parameters shown. In Fig. ^Jthis is not always 
the case but we found that there are only slight quantitative differences between the results 
of the figure and the corresponding ones with APBC The important feature to note is that 
the effect of J is to increase the peak of S(ir, it). By comparing the height of the latter peak 
for the two densities it is clear that AF correlations increase as the system approaches to 
half-filling. For n = 0.75, S(k) shows a wide structure for the usual Hubbard model, and for 
large U the peak is placed at incommensurate positions k 7^ (tt, it) (cf. lower panels of Fig. 
Hip . Remarkably, the effect of J is to shift these peaks to the AF vectors. All these features 
are consistent with the idea that J drives an enhancement of AF correlations relative to the 
usual Hubbard case. The small size of the cluster does not allows us to have an estimate of 
the AF correlation length. It can be, however, noted that for the lower density (cf. Fig. [TTj) 
S(k) spreads out on a wide range of k vectors surrounding (it, tt), consistent with a picture 
of short-range AF correlations. Instead, for lower doping, the structure evolves to a sharper 
peak around (tt, tt) suggesting larger coherence lengths. It is also interesting to note that 
the increment of AF correlations (from n = 0.75 to n = 0.875) is accompanied by a decrease 
of the pairing correlations (see Fig. 6), in agreement with the RG results from the previous 
section. 



SUMMARY AND CONCLUSIONS 



The main result of the present study is that the cooperative effect between the nearest 
neighbor exchange interaction J and the on-site Coulomb repulsion U increases in a larger 
than simply additive way the antiferromagnetic and d x 2_ y 2 superconducting tendencies in 
2D. For our analysis we used the angle resolved renormalization group including self-energy 
corrections and the exact diagonalization methods. 

We have first considered the repulsive and the attractive Hubbard models and we have 
calculated the self-energy- dressed dominant correlation functions at half-filling and at finite 
doping. In the repulsive (U > 0) case self-energy effects reduce radically all two-particle 
correlations and destroy their divergences near the critical scale. At half-filling the SC 
correlation function is below its U = value, while the AF one remains stronger than 
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its value for U = 0, but loses the divergence. This behavior, discussed in Ref. |4l|, is a 
signature of the Mott localization tendencies with simultaneous build-up of short ranged 
AF correlations. At finite doping, surprisingly and contrary to previous predictions made 
by the RG theory without self-energy corrections (Ref. 0]), the superconducting instabilities 
are strongly reduced by self-energy corrections even deeply in the BCS regime, i.e. when 
the Fermi surface is badly nested and umklapps are irrelevant. In the case of the attractive 
interaction (U < 0) our RG results are in complete agreement with previous studies 43, 44]. 
At half-filling the s-wave superconductivity and the charge density wave correlations are 
degenerate. The flow of the two correlation functions looks similar to the one of the AF 
correlation function for repulsive Hubbard model at half-filling. Just as in the repulsive 
case, the self-energy effect regularizes the flow of correlation functions near A c at half-filling 
and no phase transition occurs at this scale. The effective action of the regime below A c 
was discussed by Schulz . Contrary to the half-filled case, at finite doping the attractive 
Hubbard model shows a convincing on-set of the superconductivity. In comparison to the 
repulsive case where d x 2_ y 2 SC susceptibility is only weakly enhanced, in the attractive 
case the s-wave SC correlations at A c are not destroyed by the self-energy, while CDW 
susceptibility remains stuck to its U = value. At this point some general remarks are 
in order. We have seen that there are fundamental differences between RG flows of the 
repulsive and of the attractive Hubbard models. The fluctuations in the repulsive model 
are much stronger and probably fatal for superconductivity. They tend not only to decrease 
the magnitude of the SC susceptibility but also to decrease the energy A* = 8texp(— I*) at 
which it begins to depart from the behavior of the noninteracting one. Our present study 
is unable to say if the superconductivity is stabilized or not at some energy lower than A c . 
However, the absence of all divergences indicates that the scenarios with pre-formed pairs 
of AF and SC type are relevant even in the weak coupling limit. These RG results agree 
with the ED analysis for the intermediate-to-strong couplings, where no superconductivity 
was detected. The situation is fundamentally different in the attractive Hubbard model, 
where the superconducting instability is robust upon self-energy corrections. This is the 
RG version of the well known fact that the incoherent preformed pairs can live only at 
intermediate-to-strong interaction [46]. 

The study of the two Hubbard models was a necessary introduction to the RG analysis 
of the t — J — U model. The question that we answer is if the superconductivity of the 
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t — J — U model is enhanced or reduced with respect to the simple cases J = and {7 = 0. 
The answer is that J and U cooperate to a) increase the critical cutoff b) keeping the 
renornalized susceptibility enhanced relative to the noninteracting one. We also found that 
the latter effect is observed not only on the superconducting side of the phase diagram but 
also on the antiferromagnetic one. We thus have an evidence for the U-J synergy effect. The 
phase diagram on Fig. shows in particular the cases {U = 1.6t, J = 0), {U = 0, J = 1M) 
and (U = 1.6t, J = 1.6i). The critical temperature of the third case is much higher than 
for the first two, while the onset of d x 2_ y 2 SC remains convincingly large in the flow of the 
correlation function. 

Exact-diagonalization results support the above picture. In fact, for large enough J and 
U, we find below half-filling a clear enhancement of the SC correlations with d x 2_ y 2 symmetry 
which is accompanied by a change in the type of boundary conditions leading to the lowest 
energy. This behavior is consistent with the idea that the large Coulomb repulsion spreads 
out the Fermi surface towards sectors of the Brillouin zone where the interaction J has the 
largest amplitude in the BCS channel with d x 2_ y 2 symmetry. In such a way, we can think 
that particles are pushed to a region of the phase space, where the attractive interaction is 
most efficient to organize them into pairs. Similar arguments can be proposed to explain the 
enhancement of AF correlations at half-filling since the J interaction has components along 
AF and SC channels. In particular, in Ref. [41] it has been found that the angle-resolved 
quasiparticle weight Z{6) renormalizes in such a way that it displays a maximum in regions 
of the Fermi surface that are separated by the magnetic vector Q = (vr,7r). We have not 
found strong differences between the behavior of Z(6) in the usual Hubbard model and the 
t — J—U one neither at half-filling nor for finite doping. Therefore, the most likely scenario at 
half-filling is that the larger population of convenient regions of the k-space become available 
to be exploited by the component of the J interaction along the AF channel. 

We have not found any indication of bond-order-like instability. This is in disagreement 
with some mean-field predictions but in full agreement with other RG studies on van- 
Hove fillings of the t — t' — J — U models 0|. Our conclusion regarding the combined J — U 
mechanism to drive large superconducting correlations is also in agreement with previous 
investigations based on Fermi liquid arguments and quantum Monte Carlo simulations j^]. 
Within the repulsive Hubbard model, we obtained by renormalization group that the self- 
energy corrections are fatal for superconductivity and by exact diagonalization that the 
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superconductivity is unlikely. This result is important because it reconciles the N-patch RG 
with exact diagonalization and with a number of other approaches. Our results show that 
while the d x i_ y i pairing exists in the Hubbard model (because vertex indeed diverges), the 
onset of macroscopic superconductivity is suppressed. This result is in agreement with the 
very recent findings by Plekhanov and co-workers |47| . 

All above remarks indicate that the J interaction is vital for the superconductivity while 
the U interaction increases considerably the tendency towards pair formation. They also 
suggest that minimal microscopic models supporting recent phenomenological proposals for 

n 

colossal effects in the phase diagram of the high T c compounds [1] may be based on these 
two interactions. In fact, our results indicate that the interaction J provides a kick to the 
potential or weak antiferromagnetic and superconducting tendencies of the Hubbard model, 
that triggers a gigantic response in the system. This picture resembles the behavior of man- 
ganites where colossal magnetoresistance effects take place in response to external magnetic 
fields and this analogy is behind the proposal of Ref. The appropriate microscopic 

approach should be based on a recently developed N-patch renormalization group theory for 
disordered and interacting imperfectly nested system Qj. On the other hand, as discussed 
in our introductory section, derivations starting from the three band model for the cuprates 
also support the idea that the t — J — U model is a good candidate for a minimal one-band 
Hamiltonian for these materials. 
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